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This paper presents a model with ellipsoidal scatterers for applications to polarimetric remote 
sensing of anisotropic layered media at microwave frequencies. The physical configuration includes 
an isotropic layer covering an anisotropic layer above a homogeneous half space. The isotropic layer 
consists of randomly oriented spheroids. The anisotropic layer contains ellipsoidal scatterers with a 
preferential vertical alignment and random azimuthal orientations. Effective permittivities of the 
scattering media are calculated with the strong ftucutation theory extended to account lor the 
nonspherical shapes and the scatterer orientation distributions. On the basis of the analytic wave 
theory, dyadic Green’s functions for layered media are used to derive polarimetric backscattering 
coefficients under the distorted Bom approximation. The ellipsoidal shape of the scatterers gives rise 
to nonzero cross-polarized returns from the untilted anisotropic medium in the first-order approxi- 
mation. Effects of rough interfaces are estimated by an incoherent addition method. Theoretical 
results and experimental data are matched at 9 GHz for thick first-year sea ice with a bare surface 
and with a snow cover at Point Barrow, Alaska. The model is then used to study the sensitivity of 
polarimetric backscattering coefficients with respect to correlation lengths representing the geometry 
of brine inclusions. Polarimetric signatures of bare and snow-covered sea ice are also simulated 
based on the model to investigate effects of different scattering mechanisms. 


1. INTRODUCTION 

Technology in remote sensing has been advanced 
considerably, especially for airborne and space- 
borne radars with multifrequency and multipolar- 
ization capabilities. Theoretical models have been 
developed to interpret multifrequency data for re- 
mote sensing of geophysical media. Several ap- 
proaches including radiative transfer, modified ra- 
diative transfer, and analytic wave theories have 
been considered in model developments [Tsang et 
al ., 1985]. The analytic wave theory, which remains 
mathematically tractable for some complexity of 
the media, preserves phase information which is 
appropriate for the calculation of polarimetric scat- 
tering coefficients [Nghiem et al ., 1990]. Under the 
Born approximation for sparse and tenuous media, 
conventional scattering coefficients have been de- 
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rived for isotropic multilayered [ Zuniga et al., 
1979], anisotropic layered [Lee and Kong , 1985], 
and isotropic-anisotropic layered [Borgeaud et al . , 
1986] configurations. Fully polarimetric scattering 
coefficients have also been calculated for isotropic 
[Borgeaud et al., 1987] and anisotropic [Borgeaud 
et al., 1989] layered configurations. In these mod- 
els, scatterers are spherical or spheroidal and cross- 
polarized returns come from second- (or higher-) 
order [Zuniga et al., 1980] or tilted anisotropy [Lee 
and Kong , 1985]. For denser or less tenuous media 
the Bom approximation needs to be improved. The 
distorted Bom approximation has been applied to 
random medium modeling [de Wolf, 1971; Rosen - 
baum and Bowles, 1974; Lang, 1981]. This approx- 
imation considers dissipation loss, scattering loss, 
and the modification of wave speed due to scatter- 
ers; thus multiple scattering has been included to 
some extent. For media with strong permittivity 
fluctuations, the strong fluctuation theory is used to 
derive anisotropic effective permittivities [Tsang 
and Kong, 1981a; Stogryn , 1987]. Tsang and Kong 
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Fig. 1. Physical configuration of layered media. 



Fig. 2. Illustration of a sea ice horizontal thin section; the 
arrows represent random horizontal orientations of crystallo- 
graphic c axes and ellipsoidal brine inclusions. 


[1981a] further use the derived effective permittiv- 
ities to calculate copolarized backscattering coeffi- 
cients under the distorted Bom approximation. 
With this method, conventional backscattering co- 
efficients have been computed for an isotropic half- 
space [Tsang et aL, 1982], an anisotropic half-space 
[Tsang and Kong , 1981/?], and a configuration for an 
isotropic layer of spherical scatterers above an 
anisotropic layer of aligned spheroids [Lin et aL , 
1987]. 

Ellipsoidal scatterers have been observed in sea 
ice in the form of brine inclusions which are ori- 
ented preferentially in the vertical direction and 
randomly in azimuth [Weeks and Ackley , 1982; Gow 
et aL, 1987]. This orientation distribution effec- 
tively renders the medium anisotropic. The model 
in this paper accounts for the ellipsoidal shape, the 
medium anisotropy in a layered configuration, 
strong permittivity fluctuations, and fully polarimet- 
ric scattering of layered media. The physical con- 
figuration is illustrated in Figure 1 where % * s the 
incident angle. The covering isotropic layer is com- 
posed of randomly oriented spheroids, which have 
been used to describe scatterers in an inhomoge- 
neous layer [Nghiem et aL, 1993]. In the anisotropic 
layer, scatterers are modeled with an ellipsoidal 
correlation function with orientations characterized 
by a probability density function of the Eulerian 
rotation angles. The strong permittivity fluctuation 
theory is extended to calculate effective permittiv- 
ities, and the distorted Born approximation is ap- 
plied to derive polarimetric scattering coefficients. 


Theoretical results are compared with measured 
data at 9 GHz for bare and snow-covered sea ice. 
Effects of rough interfaces are also estimated by the 
incoherent addition approach [Lee and Kong , 
1985]. Variations of polarimetric backscattering co- 
efficients are studied for various correlation lengths 
and polarization signatures of sea ice are simulated 
and discussed for different scattering mechanisms. 

2. EFFECTIVE PERMITTIVITY 

In the isotropic scattering medium of region 1 in 
Figure 1 , embedded scatterers are modeled as sphe- 
roids, and the effective permittivity has been de- 
rived under the strong fluctuation theory [Nghiem 
et aL , 1993; Yueh et aL , 1990]. In the anisotropic 
medium (region 2), ellipsoidal scatterers have a 
preferential alignment resulting in the effective an- 
isotropy. The effective permittivity tensor of the 
scattering anisotropic medium will be calculated in 
this section for the ellipsoids with vertical alignment 
and random azimuthal orientations. 

In an inhomogeneous medium such as sea ice, sea 
water is trapped in an ice medium in the form of 
brine inclusions which are ellipsoidal as observed in 
horizontal thin sections and vertical micrographs of 
sea ice [Weeks and Ackley, 1982; Gow et aL , 1987]. 
The ice tends to grow vertically downward render- 
ing the ellipsoidal inclusions aligned preferentially 
in the vertical direction. In the absence of sea 
currents, crystallographic c axes are randomly ori- 
ented in azimuth as seen in Figure 2, which depicts 
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Fig. 3. Geometry of an ellipsoidal scattered 


a horizontal thin section of sea ice. Consequently, 
the minor axes of the ellipsoids have random orien- 
tations parallel to the horizontal plan. In this model 
a correlation function corresponding locally to a 
scatterer is used in the derivation of the effective 
permittivity tensor with the strong permittivity fluc- 
tuation theory, extended to account for the orien- 
tations of the ellipsoidal scatterers. When the aver- 
aging process is performed over orientation angle 
shown in Figure 3, the effective permittivity is an 
untilted uniaxial tensor with a vertical optic axis, 
which effectively manifests the azimuthal symmetry 
of the inhomogeneous medium. 

Let efj be the permittivity of the host medium and 
e s be the permittivity of the embedded ellipsoidal 
scatterers occupying a total fractional volume of f s . 
The subscript for the anisotropic medium (region 2) 
is omitted in this section for convenience since the 
following derivation is for region 2 only. Auxiliary 
permittivity e g = diag [e gp , e gp , e gz ] is introduced 
into the wave equation [Tsang and Kong , 1981a] for 
later consideration of the singularity in the aniso- 
tropic dyadic Green’s function G g (r 9 r). The singu- 
larity is accounted for by decomposing G g (r \ r) = 
G g (r - r) into a principal value part and a Dirac 
delta part with dyadic coefficient S which is a 
diagonal tensor with three distinctive diagonal ele- 
ments S x S y s and S z > in the local coordinates (x\ 
/, z ')• Three coordinates are related to the global 
coordinates (x, y, z) by the Eulerian rotation tensor 




■ cos 4>f 
-sin <f>f 


0 


sin 4>f O' 
cos 4>f 0 
0 1 


(1) 


Effective permittivity tensor of the inhomo- 
geneous medium is composed of a quasi-static part 
and a scattering effect part which accounts for the 
attenuation in wave amplitude and the modification 
in wave speed due to the inhomogeneities: 

e e ff=S 9 + eotf-icff-<S>r‘-!eff (2) 

where / is the unit dyad the angular brackets denote 
ensemble averaging. In (2), auxiliary permittivity e g 
and dyadic coefficient ( S ) are determined by the 
condition of secular-term elimination [Tsang and 
Kong , 1981a]. Effective dyadic scatterer £ eff , under 
the low-frequency approximation, is given by 



x 7 y,z 

d<pf p(<t>f) 2 


k,l 



dP[G g (P)] kl <t> ( (P) + [S] w 




(3) 


in which p(</>/) is the probability density function of 
orientation angle <j>f, is the variance [Nghiem 
et al., 1990], k 0 is the free-space wave number, and 
<I>^ is the Fourier transform of the normalized local 
correlation function. The anisotropic Green’s func- 
tion G g [ Tsang and Kong , 1981a], which is invariant 
under the azimuthal Eulerian rotation (1), is ex- 
pressed in the wave-vector k' domain as 


GJP) 


1 


(k’l + P;)DJP) 


lr’2 


k' x P y 


~k' x K 

Pi 

0 


(k'l + k'l)D e (P) 


U -2 


k' k' 


k'xk'y 

k'l 


TT^'i + p;n 

*9Z 


k'p 

klPeiP) 


(4a) 


where the quantities D 0 (k'), D e (k’), /c^,, and kg, are 
given by the following equations for angular fre- 
quency a) and permeability /a 0 : 


d 0 (P) = k'l + k'l + pi - k 2 gp m 
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D e (k’) = k'l + - 4 r (k'x 

Kgz 


+ k’y - k^ z ) 


Local quantities % y > , and in (7) are related to 

(4c) the elements of dyadic coefficient S = diag [Sy , S y > , 
Sy] by 


k gp = w Vo e 9P i k 2 gz = (4</) 

As in (2), the effective permittivity tensor has been 
approximated by truncating the series in the renor- 
malization method [Tsang et al ., 1985]. The validity 
condition for the approximation is |[| eff (^)]| J > w <c 1, 
where index j or m can be x, y, or z . 

The ellipsoidal scatterer is described with a nor- 
malized local correlation function of the form 




Cq 4" iSj. 1 (e 


(8a) 


= 


e — e 


9P 


e 0 *4 Sy' ( £ £ gp) 


m 




Cq 4* (c £^ 7 ) 


(8c) 


i?^(f") = exp 



where e can take on the value of e s in a scatterer or 
(5a) e b in the background medium. From (7) and (8), e gp 
and e gz can be rewritten as 


with correlation length l x > , ly, and ly in the local 
coordinates corresponding to the minor, the merid- 
ian, and the major axes of the scatterer, respec- 
tively. In this model, the correlation lengths are 
related to the effective size and shape of the scat- 
tered. This local correlation function can be re- 
duced to a spheroidal shape for two equal correla- 
tion lengths and to a spherical shape for three 
identical correlation lengths. Applying the Fourier 
transform in (5a) yields 


*t(P) = 


Ix'ly'h 


7T 2 (1 E k'lll E k']l\ E k'}l 2 z .) 2 


(5b) 


Egp Eb ^ fs 
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Eq E Egp) Cq E Sy^Efy Egp ) 

Cq E S x ‘(e s Egp) €q E S y '(€ s Egp ) 


(9a) 


Egz E b E f s 


E s Eg Z Cp 4" S Z ’(Eb Eg Z ) 

1 _ fs 4 " S Z '(E S ~ Eg Z ) 


( 9 b) 


respectively. The average dyadic coefficient (S) in 
the global coordinates is obtained by the averaging 
integration over the probability density function of 
orientation: 


For random horizontal orientations with no prefer- 
ence in azimuthal direction, the probability density 
function of orientation is simply 

p(<f>f) = 1/ (2ir) (6) 

To calculate the effective permittivity tensor ac- 
cording to (2), e g and ( S } need to be determined. 
Because of the global azimuthal symmetry, auxil- 
iary permittivity e g in coordinates (x, y, z) is 
uniaxial, as indicated previously, and the elements 
in Eg are subject to the condition 0 = 0 such that 
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( 10 ) 


Coefficients S x > , Sy , and S z > are derived from the 
secular elimination condition [ Tsang and Kong , 
1981a] which gives 



E 0 y 2 cos 2 <f> 
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Sz' = 
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In (11) the integrations over </> can be carried out 
numerically, and quantities a , y x , and y y are defined 
as 
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respectively. As observed from (14), expressions 
for /y, I y , and I z > are necessary to complete the 
derivation of the anisotropic effective permittivity. 
The result for I z > is 



-Co 


(/; + 4 ) 


(16a) 


Also, because of the azimuthal symmetry of the 
scattering medium, the effective scatterer tensor 
has the uniaxial form 


ay 2 \a y 2 
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which is obtained by substituting the Green’s func- 
tion, the correlation function, and the probability 
density function of orientation into (3). The integra- 
tions are carried out with a procedure involving 
branch cuts of complex variable functions [Nghiem 
et aL , 1990]. The following equations determine the 
diagonal elements of £ eff : 

1 

feffp = ” + + 8fy’(I y ’ + Sy)] (14 a) 


£eff 2 - 5 ^'(fz' + S z .) 


(14 b) 


where l xy l y and Iy are obtained by comparing 
(14) to (3). In (14), variance 8^, and && are 
given by 
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For Iy the result is expressed as follows: 
C2tt kl 
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The result for I y ■ can also be written in terms of /$ 
and /£ as 

k ^ 

/, = | d* - dht cos 2 1% + / 2 . r ; sin 2 4, 1 J) (18) 

Jo 77 


where the relation Iy(l x ', ly , / z ) = I x '(ly 9 l x > 9 l z ) 
has been used. Values of square-root and inverse- 
tangent functions in the above equations have been 
chosen in the respective principal Riemann sheets 
[Nghiem et al. , 1990] with branch cuts on the 
negative real axis of the complex plane. 

Substituting e g , ( S ), and £ eff in (2) yields the 
uniaxial effective permittivity tensor e eff = diag 
E^eff p , e e ffps e Qffz ] whose elements are 

£ eff p ^ gp eff p l ( 1 ^ p £ eff p ) ( 1 9fl) 


medium of permittivity e b2 . Fractional volume of 
scatterers is f sn for n = 1, 2 in region 1 and 2, 
respectively. Region 3 is the underlying homoge- 
neous half-space of permittivity e 3 . In this section 
the distorted Born approximation with permittivi- 
ties obtained in the last section is applied to derive 
the complete set of polarimetric backscattering co- 
efficients constituting the covariance matrix or the 
Mueller matrix, which characterizes polarimetric 
scattering properties of the layered media. 

Polarimetric backscattering coefficients have 
been defined with ensemble averages of scattered 
fields [Nghiem et al., 1990]. The averages are 
calculated with spatial integrations over products of 
elements of DGF G 0n , field vector F n , and fourth- 
rank correlation tensor Cg n (for n — 1, 2) as follows: 


x.y,z 
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ij,k,l,m Jo JV 2 


E eff z ^ gz ^o£effz/(l ^z^effz) (19£?) 

As seen from the above expressions, effective 
permittivity tensor e eff is anisotropic with the optic 
axis in the vertical direction. In the next section the 
anisotropic effective permittivity is used in the 
derivation of polarimetric backscattering coeffi- 
cients with dyadic Green’s functions (DGF) for the 
layered configuration. 

3. SCATTERING COEFFICIENTS 

Consider the configuration in Figure 1 for the 
following definitions of physical parameters in the 
layers. Region 0 is the upper half-space with per- 
mittivity £q. Region 1 of thickness d x is an isotropic 
medium consisted of randomly oriented spheroidal 
scatterers of permittivity e sX in a background me- 
dium of permittivity e bx . Region 2 of thickness d 2 is 
an anisotropic medium composed of ellipsoidal 
scatterers of permittivity e s2 embedded in a host 


• I dfjC £2jUm(r2 9 <t>f) 

Jv 2 

•[<Wf, ^)XF 2 ,(r 2 )>]-[(Go 2 //(r- r?)><F 2w (r?)>]* 

( 20 ) 

where V x is the volume occupied by region 1 and V 2 
is the volume occupied by region 2. The DGF and 
the mean fields have been obtained [Nghiem et al 
1990]; the correlation functions need to be specified 
next. The integrations are then carried out to derive 
the scattering coefficients. 

For region 1 , elements of correlation tensor C ^ in 
(20) have been defined for the spheroids [Nghiem et 
al., 1993]. In region 2, elements of correlation 
tensor in the global spatial domain are 

C jklm (^2 » ?2> 4>f) = (€ 2 jk (^2 )£ 2 Im 0? )| )) (21) 

To facilitate the integration of (20), Fourier trans- 
forms of correlation functions are introduced for the 
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statistically homogeneous scattering medium under 
consideration: & 2 xy 


S bl £ g2p 

e 0 4 Slx'i^bl ~~ e g2p) 


C^2jklm(^2> ^2' <t>f) ~ 
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which is expressed in the global coordinate system 
(Jc, y, z) and related by Eulerian rotation transfor- 
mation T, which is determined by (1), to the follow- 
ing nonzero correlations in the local coordinate 
system (*', y\ z ! ): 
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where Q^ip') is given in (5b), which is the Fourier 
transform of (5a), and variance values of 8 2 in (23) 
are 
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On the basis of the invariant property of the Fourier 
transform under the rotation transformation, spec- 
tral density &2jkim(P) * n the global coordinates can 
be related to those given in (24) by 

(24b) 

^2 (P) = ®&(P'x = Px cos 4> f + p y sin <f> f , 

Py = -p x sin c p f + p y cos </>/, p' z = p z ) (25) 

In the global coordinates, the rotation transforma- 
(24c) tion together with the above invariant property cast 
the anisotropic spectral densities into the form 
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jktmi.fi) - (26) 

where variance &2jklm ls dependent on the Eulerian 
angle <f>f as 
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The Eulerian rotation has been applied to arrive at 
(27) where T indicates an element of the following 
transforming tensors: 
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sin <t>f cos <t>f 
0 


sin <j>f cos <f>f O' 
sin 2 <f>f 0 

0 0 


(28a) 


derived [Nghiem et aL , 1990]. For I%${ m , the result 
is 


_ . 252 ^ 2 * (2 

*2 jklm 1 ,4 

™2 Z* 

exp [~i(K pq - K rs )d 2 ] 

i^pq ~ ^2) (^ pq ~ ^2) (* pq K rs ) 

CXp D*(^r 5 K pq )d 1] 

(^rj ^2) (^ rs ^2) ^ pq) 


~ F 2 (k 2 ) “ 02(«2) 


(30a) 


where ^ 2 (^ 2 ) an d £? 2 ( K *)* resulting from complex 
integrations, are determined by 


sin 2 <t>f —sin <f > / cos <£/ O' 
-sin (pf cos <pf cos 2 </>/ 0 

0 0 0 


(286) 


T zl =l-T xx - T yy = diag [0, 0, 1] 


(28c) 


By substituting the above correlation functions, 
the dyadic Green’s function of the layered media, 
and the mean fields into (20), the correlations of the 
scattered field can now be found. Rearranging all 
the coefficients from the DGF and the mean fields, 
the scattering coefficients can be written as 


& ILTVK 



d<t>/ <!>/) 


Plixi) 

_ i(d 2 ~ d ,) exp [-i(k 2 - k„H,] exp [<(k 2 - K pq )d 2 ] 
( 2 ilm k 2 ) 2 (k 2 _ x m ){k 2 - K rs ) 

exp [-i(K p , - K„)d,] + exp [-/(k m - K n )d 2 ] 
(2/Im k 2 ) 2 (k 2 - K pq )(K 2 - K rs ) 


exp [-i(k 2 - K rs )dj ] exp [i(k 2 - K pq )d 2 ] 
(2ilm k 2 ) 2 (k 2 - K pq )(K 2 - K rs ) 

I 1 1 

1 1 

llm K 2 K 2 r Xpq *2 - Krs 


(306) 


62 (^ 2 ) 


-1,1 x,y,z 

• E 2 KrjM?*jrJw d m + < 

a,b,c y d j,kj,m 

ou,od 

/*2 eu,ed 

■ * dd> fP2 (4> f ) 2 E (29) 

-'O p t q t r y s j,k y l,m 

The scattering contribution of the first term in (29) 
from region 1 due to the spheroids has been ob- 
tained [Nghiem et aL, 1993], and the solution con- 
sists of 16 terms for downgoing and upgoing inci- 
dent and scattered waves. Since region 2 is 
effectively anisotropic, incident and scattered 
waves can be ordinary or extraordinary and going 
down or up; these give rise to 216 terms in the 
second term of (29). All coefficients have been 


i(d 2 ~ d x ) exp [i{k * 2 - K pq )d{\ exp [-/(k* “ *rs)di] 
( 2 dm k* 2 ) 2 (k * 2 - K pq )(K 2 - K rs ) 


+ 


exp [i(k * 2 - K pq )d\] exp [-/(** - K r 5 )d 2 ] 
( 2 /Im k*) 2 (k* - k pq )(k 2 - Krs ) 


1 


dm k * 2 k * 


* 

K 2 


(30c) 


In (30) the pole k 2 is computed with the following 
equations: 


K 2 = d 2 ^\ + k\j\ x . + kl/ 2y . (31 a) 

k<t, x = k pi cos (4>j - 4> f ) - k ps cos (<t> s - <t> f ) (31 b) 


k(py k p i sin i4>! d* f) k ps sin 


(31c) 
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Corr. lengths l x ^l y * ,t z » 
Fractional vol f 3 = 4.5% 

Underlying medium : 
e = (45.0 4- i 40.00) c 0 


z — 0 m 

SEA ICE 

z — ■ 1 .65 m 
SEA WATER 


Fig. 4. Physical configuration of bare thick first-year sea ice 
with smooth surface. 


where </>; and <f> s are azimuthal angles of incident 
and scattered waves, respectively. The remaining 
integrations over Eulerian angles in (29) can be 
carried out numerically. When the low-frequency 
condition is valid, the integrations can be done 
analytically. It should be noted that (29) expresses 
the scattering coefficients in the scattered basis 
which can be transformed to the incident basis by 
changing the sign (r hTZIK and cr VThK [Nghiem et al . , 
1990]. Scattering effects of rough boundaries at 
medium interfaces can also be estimated by inco- 
herently adding the contribution from the rough- 
surface scattering with consideration of the propa- 
gation loss in the calculations of scattering 
coefficients. Compared with the model for the case 
of vertically oriented spheroids [Nghiem et al., 
1990], the new model can give higher copolarized 
backscattering coefficients due to higher total cross 
section of the ellipsoids for the same fractional 
volume. Another difference is that the new model 
provides nonzero cross-polarized backscattering in 
the first-order distorted Born term, whereas the 
former can only account for the cross-polarized 
return in the second or higher-order term when the 
spheroids are vertically oriented. 

4. RESULTS AND DISCUSSION 
4.1. Data comparisons 

In this subsection, theoretical results are com- 
pared with experimental measurements for back- 


scattering coefficient at 9 GHz as a function of 
incident angles. Radar backscatter data were col- 
lected for thick first-year sea ice near Point Barrow, 
Alaska [Kim et al ., 1984]. For the sea ice with a 
bare surface, the configuration is illustrated in Fig- 
ure 4, representing a sea-ice layer over seawater. 
The ice layer was 1.65 m thick and contained brine 
inclusions. From sea-ice characterization data 
taken at Point Barrow, Lin 11988] obtained brine 
permittivity e v1 = (38.0 4- 41.0)e 0 (refer to Vant et 
al. [1978] and Stogryn and Desargant [1985]) and 
fractional volume f s] = 4.5% (from equations by 
Frankenstein and Gardner [1967]). The background 
ice permittivity was estimated to be e b \ = (3.15 + 
/0.002)e 0 and the permittivity of seawater e 2 = (45.0 
+ /40.0)e 0 . With correlation length l lx f = 0.70 mm, 
l i y ’ = 0.25 mm, and /jv = 1.20 mm, theoretical 
results match well with the experimental data, as 
shown in Figure 5, for copolarized backscattering 
coefficient cr hh and cr vv at large incident angles, and 
cross-polarized backscattering coefficient a bv over 
the range of incident angles. 

The disagreement in the copolarized backscatter- 
ing coefficients at small incident angles is due to the 
scattering from rough surface. In the configuration 
of Figure 4, the interfaces have been assumed to be 
smooth and the rough-surface scattering has been 
ignored. Effects of the roughness at the lower 
interface of the thick first-year sea ice are negligible 
since electromagnetic waves at 9 GHz are attenu- 
ated before the seawater is reached. To account for 
rough-surface effects at the top interface, a Gaus- 
sian roughness is now considered with height stan- 
dard deviation a s = 0.6 mm and surface correlation 



Fig. 5. Backscattering coefficients from bare thick first-year sea 
ice with smooth surface. 
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length l s = 1.5 cm as depicted in Figure 6, where all 
other physical parameters are the same as in Figure 
4. The rough-surface contribution, calculated with 
the small perturbation method (SPM) [ Tsang et al. y 
1985], is incoherently added to the volume scatter- 
ing. The comparison between theoretical and exper- 
imental results for the copolarized returns is im- 
proved at the low incident angles as seen in Figure 7. 
For this surface the additional contribution to the 
cross-polarized return is small and is actually ignored 
in the SPM model applied in this case. The cross- 
polarized return therefore remains unchanged. 

The configuration for snow-covered sea ice is 
shown in Figure 8. The snow is 10 cm thick accord- 
ing to the average of the reported thickness range 
[Kim et al,, 1984]. The fractional volume of ice 
grains in snow is 34% [Lin, 1988]. In the snow layer 
the oblate spheroidal shape is assumed for the ice 
grains with correlation length l p = l z /6 = 0.15 mm. 
Physical parameters for the sea-ice layer and the 
seawater are kept unchanged. Both the top and the 
middle interfaces are rough with standard deviation 
(t s = 0.1 cm and correlation length l s = 1.5 cm. 
Theoretical results compare well with measure- 
ments from the snow-covered sea ice as seen in 
Figure 9 for conventional backscattering coeffi- 
cients. Theoretical cross-polarized returns are, 
however, lower than the measurements. For a me- 
dium with high fractional volume such as snow, a 
contribution to a hv can come from higher-order 
scattering which is ignored in the present model. 



c b — (3.15 + i 0.002) 6q 

<r 5 = (38.0 + »41.00) e 0 SEA ICE 

Corr. lengths £ x ',( y f, l z > 

Fractional vol f s = 4.5% 
z = —1.65 m 


Underlying medium : 
c = (45.0 + i 40.00) c 0 


SEA WATER 



Fig. 7. Backscattering coefficients from bare thick first-year sea 
ice with rough surface. 


Regarding the scattering mechanisms, volume scat- 
tering is dominant at large incident angles. Rough- 
surface effects are strong at small incident angles 
with more contribution from the snow-ice interface 
due to the larger permittivity contrast. Oscillations 
seen in the theoretical curves in Figure 9 are caused 



£ s = 1.5 cm (t 8 = 0.1 cm 
€ b = ( 1.00 + % 0 . 000 ) € 0 
e, = (3.15 + i 0.002) e 0 
Corr. lengths t p > and l z i 
Fractional vol f a = 34 % 


i B = 1.5 cm (X, = 0.1 cm 
€ b = (3.15 + i 0.002) € 0 
= (38.0 + t41.00) c 0 
Corr. lengths 
Fractional vol f a = 4.5% 


z = 0 m 


SNOW 


z = -0.10 m 


SEA ICE 


z = —1.75 m 
SEA WATER 


Underlying medium : 

£ = (45.0 + * 40.00) £ 0 
Fig. 8. Physical configuration of thick first-year sea ice with 


Fig. 6. Physical configuration of bare thick first-year sea ice 
with rough surface. 


snow cover. 
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Fig. 9. Backscattering coefficients from thick first-year sea ice 
with snow cover. 


by the boundary effect of the snow cover. Also, the 
theoretical oscillations in cr hh are stronger than 
those in a hv and those in a vv are weakest. This is 
because stronger reflections of horizontal waves, 
due to the medium boundaries, make the wave 
interference more effective. Comparing the cases 
with and without snow cover shows that the low- 
loss dry-snow layer can enhance the backscattering 
due to the introduction of more scatterers (ice 
grains) and a better impedance match between the 
air and the sea-ice layer. 

The above results show that the scattering from 
the inhomogeneities is dominant at large incident 
angles, while the contribution from the rough sur- 
face is important at small incident angles for copo- 
larized returns. The cross-polarized return, which 
occurs as a result of the nonspherical scatterers, 
comes from the first-order term under the distorted 
Born approximation in this model. In the next few 
subsections, effects of correlation lengths and po- 
larization signatures of sea ice are considered. 


4.2. Study of correlation lengths 

The geometry of brine inclusions in sea ice is not 
available from ice characterization data for the 
Point-Barrow sea ice under consideration. To study 
effects of the scatterer size and shape, backscatter- 
ing coefficients are calculated with variable correla- 
tion lengths. To isolate effects of correlation 
lengths, bare sea ice with smooth surface configu- 
ration in Figure 4 is considered to avoid the contri- 
bution from the snow layer and the rough inter- 
faces. The correlation lengths used in the data 
matching serve as a reference. The scatterer is 
expanded, contracted, thickened, thinned, elon- 
gated, and shortened by varying the correlation 
lengths as shown in Table 1, while the fractional 
volume is kept constant in all cases. 

When all correlation lengths are varied by the 
same factor, the scatterer shape remains unchanged 
since the correlation length ratios stay the same. 
For the same fractional volume, the effective num- 
ber of scatterers decreases as the scatterer is ex- 
panded and increases when the scatterer is con- 
tracted. Conventional backscattering coefficients as 
well as correlation coefficient p = cr hhl J ( 07 ^ 0^) 1/2 
[Nghiem et al 1990] between horizontal and ver- 
tical returns are calculated for the volume expan- 
sion and contraction with the same factor of 1.5 3 . 
The results are plotted in Figure 10. For and 
a vv , the scattering coefficients are reduced much 
more when the correlation volume of the scatterer 
is contracted as compared to the increases in the 
case of volume expansion. This shows a nonlinear 
relationship between volume scattering and the 
scatterer size. In this case an underestimate in the 
scatterer volume, compared to the overestimate 
with the same factor, will lead to more variations in 
copolarized backscattering coefficients. An inter- 
esting observation is the insensitivity of cross- 


TABLE 1. Correlation Lengths for Various Scatterer Shapes 


SHAPE 

/; 

‘y 

l'z 

I'Jl'y 


/'//;. 

I'J'yl'z 

Reference 

0.700 


■ffil 

2.800 

1.714 

4.800 


Expanded 

1.050 


in 

2.800 

1.714 

4.800 


Contracted 

0.467 

0.167 

0.800 

2.800 

1.714 

4.800 


Thickened 

0.467 

0.375 

1.200 

1.245 

2.570 

3.200 


Thinned 

1 .050 

0.167 

1.200 

6.287 

1.143 

7.186 


Elongated 

0.572 

0.204 

1.800 

2.800 

3.147 

8.824 


Shortened 

0.857 

0.306 

0.800 

2.800 

0.933 

2.614 
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Fig. 10. Effects of expanded and contracted scatterer volumes. 


polarized ratio e = cr^/cr^ as shown in Figure 10. 
This is because the depolarization is caused by the 
nonspherical shape of the scatterer, which is kept 
unchanged. For correlation coefficient p, the results 
show that horizontal and vertical waves become 
more correlated as the correlation volume in- 
creases. For larger volume of the scatterer, the sea 
ice medium is more lossy and the wave path is 
effectively reduced. Consequently, the anisotropic 
effect causing the decorrelation between horizontal 
and vertical waves is weakened. 

When correlation length ratios are varied while 
the volume of an individual scatterer and the total 
fractional volume are kept constant, the scatterer 
shape is deformed and the number of scatterers 
remains unchanged. For the same vertical correla- 
tion length l z >, horizontal correlation lengths are 
now modified by a factor of 1.5 to thicken and to 
thin the scatterer as indicated in Table 1. In Figure 


11, copolarized returns are rather insensitive, since 
the horizontal cross section of the scatterer is not 
varied (l x l z ' = 0.175 mm 2 ). It is seen in Figure 11 
that the cross-polarized ratio is strongly dependent 
on the scatterer shape. The depolarization effect is 
more prominent, as the scatterer is further de- 
formed from the spheroidal shape. For the magni- 
tude of p, the variation is larger at small incident 
angles where the scatterer shape is more different as 
compared to the cases at large incident angles. The 
phase of p started from 0° at normal incidence due 
to the azimuthal symmetry [Nghiem et al . , 1992] 
and increases with incident angles. Also, the abso- 
lute value of the phase is larger for the thicker 
scatterer since the anisotropy is stronger. 

While the correlation length ratio l x 'Uy, the num- 
ber of scatterers, and the fractional volume are not 
changed, the scatterer is elongated and shortened 
by varying vertical correlation length l z > by a factor 
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Fig. 1 1. Effects of thickened and thinned scatterer shapes. 


of 1.5. Figure 12 presents results for these cases 
compared with the reference case. There are large 
changes in copolarized returns due to the corre- 
sponding differences in the horizontal cross section 
of the scatterer. The variation in the cross-polarized 
ratio is small, since horizontal correlation length 
ratio l x 'lly’ is not varied. The magnitude of p differs 
more at large incident angles, where the shape 
difference is more significant. The phase of p is 
closer to zero as the scatterer becomes shorter and 
the sea-ice medium approaches the isotropic condi- 
tion. As shown in this subsection, scatterer volume 
and horizontal cross section affect the level of 
copolarized returns. The depolarization effect char- 
acterized by cross-polarized ratio e is controlled by 
the scatterer shape. The ellipsoidal shape of the 
scatterer with a preferential alignment in the verti- 
cal direction is also responsible for the anisotropy 


of the medium and, consequently, the behavior of 
complex correlation coefficient p between horizon- 
tal and vertical returns. 

4.3. Polarimetric simulation 

To investigate how polarization signatures relate 
to the corresponding covariance matrices and con- 
vey information regarding different layered config- 
urations and scattering mechanisms, polarimetric 
data are simulated for sea ice with a bare surface 
and with a snow cover at different incident angles. 
Definitions and characteristics of covariance matrix 
and polarization signature can be found in refer- 
ences [Nghiem et aL, 1990; van Zyl and Zebker , 
1990]. The configurations under consideration are in 
Figure 4 for bare sea ice with a smooth surface, in 
Figure 6 for bare ice with a rough surface, and in 
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Fig. 12. Effects of elongated and shortened scatterer shapes. 


Figure 8 for snow-covered sea ice. As a result of the 
azimuthal symmetry of the scattering configuration, 
the covariance matrix has the form of [Nghiem et 
at . , 1992] 


C = 


- 1 0 
0 e 

_p*\y 0 


pH 

o 

y 


(32) 


where the copolarized ratio is y = a vv lu hh . Note 
that covariance matrix C (without subscript) should 
be distinguished from correlation tensor C which 
is tagged with subscript £ corresponding to corre- 
lation function and subscript n for region n - 1, 
2. Elements of the covariance matrix arse used to 
calculate the Mueller matrix and to obtain the 
polarization signature o(a r , p n a h /3,-) [Nghiem et 
at 1990; van Zyl and Zebker , 1990] for polarization 
angles a r and a,* and ellipticity angles p r and j8 f of 


the received (r) and the incident (/) waves, respec- 
tively. For a copolarized signature, a r = a, = a and 
Pr = Pi = P- 

Normalized copolarized signatures denoted by a n 
[Nghiem et at., 1990] for bare sea ice with smooth 
and rough surfaces at incident angles of 20° and 45° 
are presented in Figure 13. At % = 20°, signatures 
for smooth sea ice due to volume scattering and for 
rough sea ice due to surface scattering are similar 
because both covariance matrices have y close to 
unity, small e , and p with small phase as seen in 
Figure 14. However, |p| shown in Figure 14 is 
distinctively smaller for the volume-scattering 
mechanism with smooth surface. At % = 45° for 
both smooth- and rough-surface cases, volume scat- 
tering is dominant and correlation coefficients with 
larger phase give rise to signatures in Figures 13c 
and 13d with more distortion compared to those at 
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Fig. 13. Polarization signatures of bare thick first-year sea ice with ( a ) a smooth surface with an incident 
angle (0 O ,) of 20°, (b) a rough surface at incident angle of 20°, (c) a smooth surface with an incident angle of 45°, 
and (d) a rough surface with an incident angle of 45°. 


20°. For snow-covered sea ice, Figure 14 reveals 
that |p| is closer to unity, and the phase is smaller at 
incident angles of 20° and 45°. Thus the correspond- 
ing signatures for the snow-covered sea ice in 
Figure 15 have much less distortion due to the 
rough surfaces and the isotropic snow layer which 
mask the anisotropic information from the lower 
sea-ice layer at the X band frequency. Boundary 
effects due to snow cover can also be seen with the 
oscillations in the phase of p in Figure 14. 


5. SUMMARY 

In this paper a model has been developed to 
account for ellipsoidal shapes of scatterers in aniso- 
tropic layered media. Ellipsoidal scatterers with a 
preferential alignment in the vertical direction have 
been observed in horizontal and vertical sections of 
sea ice. Effective permittivities are derived with the 
extended strong fluctuation theory, and fully pola- 
rimetric backscattering coefficients are calculated 
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Fig. 14. Correlation coefficient p of thick first-year sea ice with smooth surface, rough surface, and snow cover. 


under the distorted Born approximation. In this 
model the cross-polarized return is obtained from 
the first-order approximation due to nonspherical 
shapes of the scatterers. Effects of rough surfaces 
are estimated with the incoherent addition method. 
Theoretical results compare well with measured 
data at 9 GHz for thick first-year sea ice with a 
bare surface and a snow cover. Changes in scatter- 
ing coefficients due to the scatterer geometry are 
investigated by varying correlation lengths which 
modify the shape or the volume of the scatterer. 
This sensitivity study also provides an estimation 
of variations in polarimetric scattering coefficients 
due to uncertainties in the scatterer sizes and 
shapes. Polarization signatures are simulated for 



Fig. 15. Polarization signatures of thick first-year sea 


the thick first-year sea ice with and without snow 
cover. These signatures convey information on 
the medium anisotropy, which are influenced by 
rough-surface or snow cover effects. Volume- and 
surface-scattering mechanisms are discussed in 
terms of conventional backscattering coefficients 
as well as polarization signatures and correlation 
coefficient p between vertical and horizontal re- 
turns. Further medium complexity such as multi- 
species, size and shape distributions, and higher- 
order volume or surface scattering have not been 
incorporated in this model and can be considered 
for future model developments with comparisons 
to well-calibrated measurements of polarimetric 
covariance matrices. 



ice with snow cover at incident angles of 20° and 45°. 
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